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In this article, we show that the exact two-body problem can be replaced by quantum jumps 
between densities written as D — \9 a ) where |$ ) and l^b) are vacuum for different quasi- 
particles operators. It is shown that the stochastic process can be written as a Stochastic Time- 
Dependent Hartree-Fock Bogoliubov theory (Stochastic TDHFB) for the generalized density 1Z 
associated to D where 1Z 2 = 1Z along each stochastic trajectory. 



I. INTRODUCTION 



The goal of this article is to prove that the dynamics of 
fermions interacting through a two-body interaction can 
be transformed into a stochastic process in the Hilbcrt 
space of Hartree-Fock Bogoliubov (HFB) states. Such a 
derivation is motivated by recent studies dedicated to the 
structure of nuclei. In nuclear physics, mean-held theo- 
ries like Hartree-Fock and HFB already provide a good 
approximation of static and dynamical properties [l| . It 
also turns out that a deep understanding of nuclei re- 
quires the introduction of correlations beyond mean-field. 
Large theoretical developments are devoted for instance 
to Generator Coordinates Methods (GCM) Q. In that 
case, some collective degrees of freedom are selected and 
the correlated ground state is constructed as a super- 
position of mean-field states (HF or HFB). It has been 
shown that the description of nuclear systems is greatly 
improved if pairing correlations are already accounted 
for, i.e. if the GCM is performed with HFB many-body 
states. Such a method has been successfully applied to 
nuclear structure when only few degrees of freedom are 
selected. However, applications of GCM techniques when 
many collective degrees of freedom are important are still 
numerically intractable. 

Monte-Carlo techniques appear as an alternative way 
of treating correlation beyond mean-field. Shell model 
Monte-Carlo theoryQ is an example of such technique. 
Recently, starting from the Hartree-Fock theory, new for- 
mulations 3,Sl3 have been proposed that combines the 
advantages of Monte-Carlo methods and mean-field the- 
ories. In that case, the exact evolution of fermionic (or 
bosonic) systems is replaced by an ensemble of stochastic 
mean-field evolutions. A possible improvement of such 
theory, which might be of great interest in the nuclear 
context, is to include pairing correlations in the trial set 
of wave functions, i.e. to consider quantum jumps be- 
tween HFB states instead of HF states. A first step in 
this direction as been made in ref. where stochastic 
dynamics between BCS states were introduced. In this 
article, we show that the dynamics of fermions interact- 



ing through a general two-body Hamiltonian 

ij ijkl 

can be mapped into a quantum jump process between 
HFB states. Here, c + operators correspond to creation 
operators associated to a complete single-particle basis, 
and V12 matrix elements are antisymmetrized. 

In the following we first introduce quantities associ- 
ated to " densities" written as a dyadic of two HFB state 
vectors, i.e. D — \^f a ) (^>b\- The flexibility of stochastic 
methods allows to consider densities with specific helpful 
properties which are precised in the first part of this work. 
Then a TDHFB equation is derived for D when correla- 
tions beyond mean-field are neglected. Finally, the full 
stochastic theory that accounts for all two-body effects 
is derived. 



II. PRELIMINARY RESULTS AND NOTATIONS 

In quantum Monte-Carlo approaches starting from an 
initial density D — |\&) (*&\, the exact system evolution is 
recovered by averaging over an ensemble densities written 
as a product of two different state vectors 

£> = |*«> <*&l • (2) 

The use of two different states is at the heart of exact 
stochastic methods. The advantage of these approaches 
is that states entering in D correspond generally to a spe- 
cific class of trial wave-function. In previous applications, 
these states have been chosen as Hartree-Fock states|5j. 

Given a specific choice of trial wave-functions, it turns 
out that stochastic reformulation is generally not unique. 
This flexibility might be used for instance to optimize 
quantum jumps and reduce the number of paths (see for 
instance [8j)- Here we will consider HFB state as trial 
state vectors and use this flexibility in a different way. 
Indeed, it turns out that the two-body problem can be 
reformulated as a stochastic process imposing additional 
relations between the states \^ a ) and \^b) along each 
path. These additional constraints, given below, lead to 
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simplified equations and derivations without restricting 
the exactness of the formulation. 



A. Choice of a subclass of densities 

Let us assume that D has the form J2J) where \^ a ) 
and | \l/ j>) can be written as a product of quasi-particle 
operators, i.e. |* Q ) = Ha+ |0) and \$> b ) = Hb+ |0). D 
will be referred in the following as a density although it 
does not necessarily meet all required properties to be 
considered as a density matrix. For each quasi-particle 
operator, two sets of single-particle wave- functions, de- 
noted by |a aj b) and \dt a ,b), are introduced. They define 
the transformation between quasi-particle states and a 
complete set of particle states as 



or equivalently: 



(fife 



\i) + (i | a a ) c+, 
+ \ a b )cf. 



(3) 



Note that, we can recover the matrix notations U, V often 
used in the HFB theory through the relation a ' U{ a = 
(a a ,b I i) and^'K,, = (i \ a a . b ). 

As usual [J, we introduce vector notations a = 
{a, a + }, b = {b, b + } and c — {c, c+}. Above linear trans- 
formations can then be written as matrix transformations 
a = W^c and b = W^c. In opposite to the standard 
HFB theory, we do not impose the transformations to be 
canonical but instead restrict ourselves to a subclass of 
quasi-particles and vacuums having two specific proper- 
ties. We first assume that 



w a m + = w a + m = 1, 



(4) 



which gives the inverse transformations c = V\4a = W a b. 
As a consequence, although the a and a + operators 
(respectively the b and b + ) do not necessarily fulfill 
fermionic anti-commutation rules, because of 10} we have 



la, °I3\ + 



= &£]+ = 0, [a a ,bt]+ = S a0 . (5) 



The second important assumption is that |\l/ a ) and \^b) 
are both vacuum for all a a and b a . As we will see in the 
following such properties might occur without any simple 
relations between the two sets of annihilation operators. 

We introduce the generalized density matrix TZ a b de- 
fined as 



Ttab 



<* 6 |c+ Ci |* a ) (* 6 \CiCj\ * a ) 



From the two assumptions, it can be shown that 



W?KW a = L 




1 



(6) 



(7) 



< a a bp > = < a+bp >=< >= 0, 



< a a bj } > 



This again can be seen as a generalization of the HFB 
case and implies TZ^b ~ T^-ab- In addition, the generalized 
density lZ a b takes a simplified form compared to the one 
generally obtained for transition densities Q. Here, we 
have 



7t a b 



Pab K ab 
* 1 T 
~ K ab 1 -Pab 



(8) 



where p£ b denotes the transposed matrix of p a b- In the 
following, to simplify notations we will omit the sub- 



script 



Different operators matrix elements can be 



expressed as 

P 

1-P 



X>«>KI> ( 9 ) 

a 

X>a}v56|, (10) 

OL 

2J \a a a b ) = - ^ \®aUb) , (11) 



with the convention — ^ Q (ij \ a a 6i b ). 

Finally, we will also use the notation | a,b H / Q ) and 
|°' b K*) (taken from ref. We have in particular 

n = ^2\ a w a )( b w a \, 
l-n = £r%X 6 v a |, (12) 

a 

with ( a W p I b W a ) = ( a V p I b V a ) = 5 a0 and 
( a Va I b W a ) = 0. This completes the different proper- 
ties associated with the subclass of densities considered 
here. 



B. Expression of the Hamiltonian and generalized 
TDHBF equation 

Using the previous properties, the action of the two- 
body Hamiltonian on the vacuum \^ a ) can be recast as 



H\* a ) = {(H)+h L +Hf- es }\y a ), 



(13) 



where we have used the compact notation (H) = 
(^b \H\ fy a ) and where is a one-body effective Hamil- 
tonian given by 
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h and A correspond respectively to matrix elements 

hij = T i3 + (i\Tr 2 (v 12 p 2 )\j) , (15) 
Ajj = ^ ^2 (ij \v 12 \ kl) K kh (16) 

hi 

which will be called mean-field and pairing field in anal- 
ogy to HFB theory. Note that expression l|14|) dif- 
fers from the one generally obtained in HFB using the 
Wick theorem because of the coexistence of two sets of 
quasi-particle operators. Starting from (|14|1 - the effective 
Hamiltonian can be recast as 

h L = \(c+ c)(l-K)Wl^+y (17) 
where Ti stands for the generalized HFB Hamiltonian 

m 

K-(_*. _tr). (18) 

We will see that expression Ijl7|l is central for further de- 
velopments. The last term of equation (|13|1 . called here- 
after residual Hamiltonian reads 

H res = \ E \Vl2\5ala) a+<$b+b+. (19) 

Performing a similar decomposition of (*S?b \ H leads to 
(* b \H=(* b \{(H)+h R + H* s }, (20) 

with 

h R = \{c+ c)K'H(l-K)(° + y (21) 

while 

H res = ^ E ( a b@b\vi 2 \5 a j a ) a a af3b 7 b s . (22) 

C. Evolution of the generalized density 1Z 

Starting from the initial density (0), the evolution of 
the system is considered assuming first that the effect 
of the residual interaction can be neglected,. After one 
time-step, due to the one-body nature of h^, the state 
\^> a + d^a) = e~ hh \^ a ) is a vacuum for the new quasi- 
particles a' a = a a + da a = e^ hh a a e~^ hL . Similarly, 



I 

(^6 + d^>b\ = (^b \ e~^ ha is a vacuum for the new quasi- 
particles b' a + = b+ + db+ = eiK hR b+e-^ hR . 

Since the residual interaction is neglected, all the in- 
formation on the system is contained in the evolution of 
TZ. From standard rules of creation-annihilation operator 
transformations 0E3>E1> we obtain: 

[h L ,c] = -(l-n)HKc, (23) 

[h R ,c\ = -nn(i-n)c. (24) 

With the help of the above anti-commutation relation- 

, , dt 1 1 dt l q dt l dt 1. 

ships, we can express e i« nt ce + is 1 ande>» s ce 
and deduce the evolution of TZ. Using the fact that ini- 
tially TZ 2 = TZ and (l-K)K = 0, we end with 

in— = (i-n)wi-Kn(i-K) 

= [H,K], (25) 

which is nothing but a TDHFB equation generalized to 
densities given by eq. Without going into further 

details, it can be shown that the density TZ fulfills all 
properties listed above and thus remains in the subclass 
of densities previously described. Therefore, the consid- 
erations made for single time step can be extended to the 
long-time dynamics. 

III. INTRODUCTION OF QUANTUM 
MONTE-CARLO METHODS 

In the previous section, we have introduced general 
properties of densities given by eq. @ which will be 
helpful for the forthcoming discussion. In addition, we 
have shown that the dynamics reduces to a TDHFB-like 
equation when the residual interaction is neglected. The 
aim of this section is to show that the residual interac- 
tion can be treated by introducing stochastic processes 
between densities described in section HTl 

A. Separable residual interaction 

In the following discussion, we concentrate first on the 
evolution of \*5f a }, keeping in mind that everything can be 
transposed to (^b\- Following the Stochastic mean- field 
approach 0, H, El , we consider that the residual part of 
the interaction can be written as a sum of of separable 
interactions in the particle-hole channel: 

(ctaPa K2I Sbjb) = - E ("o \O m \5b) (Pa |Om|7b) >( 26 ) 
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where O m corresponds to a set of single-particle opera- 
tors. Using relation (|26[) . the residual interaction H^ es 
can be recast as 



(27) 



where the set of one-body operators B m are given by 

B%=J2(*a\Om\Pb)aib+. (28) 

a/3 

Guided by previous section, we write it as 

B£ = \ ( c c+ ) (1 - ttJB^ft ( C c + ) , (29) 

where we have introduced the matrix : 
O m 



-02 



(30) 



Once -ff,P es is written as l|27|) , the introduction of stochas- 
tic process is rather straightforward. Introducing a set 
of stochastic variables d£J^ (which follow Ito rules of 
stochastic calculus 12]) with mean values equal to zero 
and variances satisfying 



(31) 



Here the ( n ) exponent stands for a specific realization of 
the stochastic process. In the following, it will sometimes 
be omitted to simplify notations. The evolution of \$f a ) 
associated to H can then be written as an average over 
stochastic evolutions in the Hilbert space of HFB state 
vectors: 



= e «w e i^+£ m ^ B £> ) (32) 



= w„ n) (* + <&)). (33) 

In this equation, same conventions as in ref. 0, are 
used and ^a^(t + dt)) correspond to different vacuum 
states. Introducing the notation 

dt. 



dSp h 



ih 



E 



(34) 



,T,(™)\ 

Wo ) is a vac- 



according to the Thouless theorem, each 

mini for the quasi-particle a a = e p h a a e p h . 

One can finally note that the stochastic evolution 
of | Wo) should be completed by an equivalent stochas- 
tic evolution for (w&|. The associated propagator and 
stochastic variables are respectively denoted by dS R h and 
d£ R . These variables should be taken statistically inde- 
pendent of to properly account for the exact dynam- 
ics of D. More explicitly, we have 



ds^ h = ~h R +j2^X, 



(35) 



where 



Bt l = ^(c c+)TU3£{\-K) 



(36) 



and = -<W'§ 



IV. NATURE OF THE STOCHASTIC PROCESS 

Similarly to the previous case, where the residual in- 
teraction was neglected, we do expect that along each 
path, the stochastic evolution of the system reduces to 
the stochastic evolution of TZ. The explicit form of the 
stochastic evolution of TZ can now be obtained using the 
commutation relationship 0, Hoi Hlj : 



e dS pp , c 



e -( 1 -^)[f n+B L ]n c 



in 



TZc, (37) 



while 



e w,c 



= p-^hf n+B R ](i-n) 



= c-TZ 



in 



where B L and B R stand for 



B L/R 



(1 - TZ)c, (38) 



(39) 



Note that equations l|37|l and l|38|) are exact thanks to 
the (1 — TZ) term. Similarly, as in previous section and 
using expression (|37|l and l|38|) . one gets the stochastic 
evolution of TZ: 



dTZ 



dt 



[H,TZ] + (l-TZ)B L TZ + TZB R {l-TZ). (40) 



Such a stochastic process, called hereafter Stochastic TD- 
HFB, is similar to the Stochastic mean- field dynamics 
|l4j except that the mean-field and normal densities are 
now replaced respectively by the HFB Hamiltonian Ti 
and density TZ. Starting from TZ = £ Q \ a W a ) { b W a \, 
evolution of TZ can be replaced by the set of equations 



\d a W a ) = {§H + (l-TZ)B L }\ a W a ) 
(d b W a \ = ( b W a \{-§-H + B R (l-TZ)} 



(41) 



Above expressions show that if initially fullfilled, the 
property (""Wp | b IU Q ) = 6 a p is true all along each 
stochastic path. Again, it can be shown that all the prop- 
erties of the class of densities considered in section UTI arc 
preserved during the stochastic evolution (J3DJ. There- 
fore, we only have to initiate the quantum jump process 
with a density which satisfies the properties described in 
sectionlH] This is the case if we start from an initial HFB 
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density D — \9) (\&|, which is the most convenient in 
practice. Note finally that the explicit form of the quasi- 
particle evolution can directly be obtained from eq. lBT|l 
while the stochastic evolution of p and k can be deduced 
from IpEUl) . 

A. Alternative form and pp — hh separable 
interaction 

In the previous section, we have developed quantum 
diffusion processes between HFB states assuming expres- 
sion l|26(l . However, recent studies in nuclear structure 
fl5l | support separable interactions in the particle-particle 
and hole-hole channels. For completness, we introduce a 
set of one-body operators G m and assume that the resid- 
ual interaction now reads 

{a a p a \v 12 \5 blb ) = -^(aa|G m |/3 h )<<5 a |G m | 7 fc}*{42) 

m 

where G m should be a skew matrix (i.e. G^ n — ~G m ) to 
respect the antisymmetrization of v\2- The formulation 
of the two-body problem can equivalently be done start- 
ing from eq. lH^ . The final result is that TZ still obeys 

a stochastic equation with similar form as (|40|l . where 
qL/r nQW rea( j s 

B L/R = £ d V %«B% + i (dr,^) * fi** (43) 

m 

and drj m corresponds to stochastic variables with mean 
value zero and 



^W')* = S mm'^, (44) 

while all other second moments are equal to zero. 
and B^ are given by 
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V. CONCLUSIONS 

In this work, we have shown that the exact dynam- 
ics of interacting fermions can be replaced by a Monte- 
Carlo method in the Hilbert space of Hartree-Fock Bo- 
goliubov states. In order to prove this reformulation, we 
have used an intermediate result, considering densities 
D = l^a) (^hl with specific properties. Neglecting the 
residual interaction, the evolution of D leads to a TD- 
HFB equation for TZ. We then have proven that the in- 
troduction of correlations beyond the mean-field picture 
can be replaced by a Stochastic TDHFB equation for TZ, 
generalizing the stochastic mean-field approach It 
should be noted that the reformulation is not unique and 
the selection of a sub-class of D is not absolutely nec- 
essary. However, in that case derivations and stochastic 
equations are more complicated. 

The stochastic theory presented here is not restricted 
to dynamical problem and could also be useful to study 
static properties of interacting systems [jg. In that case, 
real time propagation is replaced by imaginary time evo- 
lution. Monte-Carlo methods has the advantage of not 
requiring an a priori knowledge of the relevant collective 
degrees of freedom and can eventually be used as an alter- 
native to GCM. It should however be noted that Monte- 
Carlo methods still require large numerical efforts. Work 
is actually in progress to combine advantages of GCM 
and Monte-Carlo techniques. 

Finally, we would like to mention that the above the- 
ory gives an indirect proof of the fact that densities de- 
scribed in section[n]form an over-complete set of densities 
to treat the two-body problem. This might be of great 
interest even for non stochastic methods which treat cor- 
relations beyond mean-field. 
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